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Aims. The aim of the paper is to understand the emission from the photon dominated regions in Cepheus B, estimate the column densities 
of neutral carbon in bulk of the gas in Cepheus B and to derive constraints on the factors which determine the abundance of neutral carbon 
, relative to CO. 

Methods. This paper presents 15'xl5' fully sampled maps of [Ci] at 492 GHz and 12 CO 4-3 observed with KOSMA at 1' resolution. The new 
OO observations have been combined with the FCRAO 12 CO 1-0, IRAM-30m 13 CO 2-1 and C ls O 1-0 data, and far-infrared continuum data from 
vO HIRES/IRAS. The KOSMA-r spherical PDR model has been used to understand the [Ci] and CO emission from the PDRs in Cepheus B and 

to explain the observed variation of the relative abundances of both C° and CO. 

Results. The emission from the PDR associated with Cepheus B is primarily at Vlsr between -14 and -11 km s '. We estimate about 23% of 

the observed [C n] emission from the molecular hotspot is due to the ionized gas in the H n region. Over bulk of the material the C° column 

density does not change significantly, (2.0 ± 1.4) 10 17 cirT 2 , although the CO column density changes by an order of magnitude. The observed 

C/CO abundance ratio varies between 0.06 and 4 in Cepheus B. We find an anti-correlation of the observed C/CO abundance ratio with the 

■ observed hydrogen column density, which holds even when all previous observations providing C/CO ratios are included. Here we show that 

Oh this observed variation of C/CO abundance with total column density can be explained only by clumpy PDRs consisting of an ensemble of 

clumps. At high H2 column densities high mass clumps, which exhibit low C/CO abundance, dominate, while at low column densities, low 

5^ ■ mass clumps with high C/CO abundance dominate. 
-4— > ■ 
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£3 , 1 . Introduction abund ances of C° and C + in the gas phase jBoeer & Sternberg! 

■ - - ' [2005). PDRs may therefore be seen as the most ubiquitous 

Photon Dominated Regions (PDRs) are predominantly neutral, phase of ^ ISM ^ observations of which together with proper 

atomic and molecular regions where the physical and chem- modeling provide insight into the physical conditions within 

ical processes are dominated by Far Ultraviolet (FUV) radi- the molecular clouds Qn a more global sca ^ PDRs provide 

ationl [cf. I Sternberg&Dalgarno l tl99j;| Horienbach & Tielens | direct insight Mq the chemical evolution of the ISM deluged 

11223) and the major carbon-bearing species changes from C+ wkh pjjy photons and PDR models suggest rad iation-induced 

through C° to CO. While FUV photons (6 eV< hv < 13.6 eV) feedback mec hanisms that may regulate star formation rates 

are primarily responsible for the heating of the cloud surfaces and ^ column denshy of gag through giant molecular clouds 
via photoelectric effect on dust grains, at larger depths of about 

1 f\ . .1 1 1 • • , ■ , . • ■ ■ . To the south of the Cepheus OB3 association of early- 

10 mag into the clouds cosmic ray induced heating dominates. F 3 

■ . , , f 1 1 1 j t type stars lies the giant Cepheus molecular cloud, at a dis- 

Owing to the clumpy nature of the molecular clouds, contrary Jr 6 1 t - L — ■■ 1 

to the initial ideas, PDRs are not strictly confined to the sur- tance of about 730 P c Blaauw || l964 . C epheus B, wifh asize 

face of the molecular clouds. Rather the extended [Cn] and of 4 P c ' is the hottest CO component jSar R enfll97l | l979| > 

r „ -. ■ ■ , , t „ 1 1 1 j and is located near the northwestern edge of the cloud, adia- 

[C 1] emission observed from allover the molecular clouds sug- 1 ^^-fi -j J 

gest formation of PDRs deep inside the molecular clouds, at cent t0 the Hn re 8 ion S155 < Sharpless| | l959|) . The interface 

e e , • , . j ... CI ■ . • between the molecular cloud and the OB stars, is delineated 
the surfaces of clumps inundated with FUV photons escaping 

c • 1U , t ,, by the optically visible Hn region S155, whose very sharp 

into the cloud. Further, the chemical abundances of several key 3 F 3 & 3 F 

species like e.g., CN and HCN depend to a large extent on the ed § es indicate the P resence of ionization fronts bound " 

ing the dust/molecular cloud. The OB association itself con- 

W offprint requests to: B. Mookerjea, e-mail: sists of two subgroups of different ages, with die youngest lying 

bhaswati@astro.umd.edu closer to the molecular cloud JSargenll979h . lFelli et alJ dl978t) 

* Present address: Department of Astronomy, University of presented radio continuum observations which confirmed that 

Maryland, College Park, MD 20742, USA the molecular cloud Cepheus B is surrounded by an ioniza- 
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tion front driven by the radiation field of the bright stars that 
are part of the youngest generation of the Cepheus B OB3 as- 
sociation. The structure of the diffuse ionized gas and its en- 
ergy balance can be accounted for by the UV radiation com- 
ing from the brightest members of the OB3 association (viz., 
HD217086 and 217061). The physical association between the 
Cepheus B molecular cloud and the S155 Hn region was also 
co nfirmed by the H?CO a n d recombination line observations 
of lPanada & ThurrJ 1l98ll) . lMinchin et alJ Jl992h used a com- 
bination of low-J CO and far-infrared (FIR) continuum obser- 
vations to show that the Cepheus B molecular cloud terminates 
at its north-west edge with a hotspot. The hotspot lies V to the 
east of the peak in CO column density. The high temperature 
of the hotspot and the associated FIR emission cannot be ex- 
plained by only consideri ng the stars in the OB 3 association. 
It was further proposed by Minch in et al.l ([1,992) and later con- 
firmed us ing radio and nea r-infrared (NIR) continuum observa- 
tions bv lTesti et al 1 11 199.4 that the hotspot is indeed energized 
by a highly reddened main sequence star (NIR- A) of spectral 
type B1-B0.5, together with a small cluster of young stars. 



In a firs t study with the KOSMA 3m telescope, 
Beut her et all koOO) mapped the Cepheus B cloud in 2-1 and 
3-2 transitions of 12 CO, 13 CO, and C 18 Q at 2' resolution and 
used PDR models of lStorzer et alJ (l2000l) to interpret the emis- 
sion at four selected positions. The hotspot emission indicated 
the presence of shocks. Based on the derived local volume den- 
sities of ~ 2 10 4 cm' 3 , and ave rage volume density less than 
10 3 crrT 3 . iBeuther et all 1^2000) concluded that Cepheus B is 
highly clumped with clumps filling only 2% to 4% by volume 
of the cloud. In this work [C 1] at 492 GHz was detected from 
two selected positions and detection of [Ci] emission from 
deep inside the cloud was also put forward as an additional 
evidence for the dumpiness of Cepheus B. 



In clumpy or homogeneous PDR scenario both [Ci] and 
12 CO 4-3 emission arise from regions closer to the surfaces 
of the FUV irradiated clumps or molecular cloud surfaces re- 
spectively. Thus the [C 1] and 12 CO 4-3 emission trace similar 
regions. While [C 1] emission necessarily implies photodisso- 
ciation of molecular CO, the higher kinetic temperature (the 
upper energy level being 56 K) required to excite the 12 CO 4- 
3 line suggest origin of this emission from the outer surfaces 
of dense (since the critical density of the 12 CO 4-3 transition 
for collisions with H2 is ~ 10 5 cm -3 ) clumps embedded in a 
UV radiation field. Thus in order to trace the PDR we have ob- 
served largescale (15' x 15') fully-sampled maps of the [Ci] 
and 12 CO 4-3 emission from the Cepheus B molecular cloud. 
At a common resolution of 1', this dataset has been combined 
with FCRAO 12 CO 1-0, IRAM-30m 13 CO 2-1 and C 18 1-0 
data, and FIR continuum data from HIRES/IRAS. This allows 
for a detailed analysis of the PDRs in Cepheus B, since ob- 
served transitions span a substantial range of critical densities 
and temperatures (Table 0. 



Table 1. Parameters of the observed transitions used to study 
the PDRs in Cepheus B. 



Species 


Transition 


E u /k 








(K) 


(cm- 3 ) 


[Ci] 


3 Pi ^ 3 Po 


23.6 


1 10 3 


12 CO 


1-0 


5.6 


4.8 10 3 


13 co 


2-1 


16.8 


3.2 10 4 


12 co 


4-3 


56 


2.6 10 s 



2. Datasets 

2.1. KOSMA Observations of PDR tracers: [C\] and 
12 CO 4-3 

We have used the KOSMA 3m submillim eter te l escope 
on Gornergrat, Switzerland JWinnewisser et all ^986; 
iRramer et alJll998i) to observe the emission of fine- structure 
line of neutral carbon at 492 GHz (609 /urn, 3 Pi- 3 Po; hereafter 
[Ci]) and J = 4-3 transition of 12 CO emission at 461 GHz 
from the PDRs associated with Cepheus B. We used the 
Submillimeter Array Receiver for Two frequencies (SMART) 
JGraf etal. l2002t> on KOSMA for these observations. SMART 
is a dual-frequency eight-pixel SIS-heterodyne receiver ca- 
pable of observing simultaneously in the two atmospheric 
windows centered around 650 yum and 350 fim. The IF signals 
were analyzed with array-acoust o-optical spectro meters with a 
spectral resolution of 1.5 MHz faorn et all 19 99). 

The observations were performed in the position-switched 
On-The-Fly mode between November 2004 and March 2005. 
The HPBW of the [Ci] and 12 CO 4-3 observations are 60" 
with a beam efficiency of r] m \, = 50%. Atmospheric calibra- 
tion was done by measuring the atmospheri c emis sion at the 
OFF-position to derive the opacity (Hivama 1998). Sideband 
imbalances were c orrected using standard atmospheric models 
JCernicharo 1 1 985i) . All data presented in this paper are in units 
of main beam temperature (T m b), calculated from the observed 
calibrated antenna temperature (TV) using the derived beam ef- 
ficiencies (?7 m b) for the relevant telescope, T m b = T* /r] m ^,. From 
observations of reference sources we estimate the relative cal- 
ibration accuracy to be ~ 15%. Data was reduced using the 
GILDAS 1 astronomical data reduction package. 

2.2. Further CO data 

We have complemented the [Ci] and 12 CO 4-3 datasets for 
Cepheus B with low-J CO observations available in literature. 
We have obtain ed the 12 CO 1-0 da taset from the FCRAO outer 
Galaxy survey faever et all 19981) and 13 CO 2-1 and C 18 Q 1- 
data sets observed with the IRAM 30m by llJngerechts et alJ 
(2000) in 1995. The beam efficiency, rj m b (defined as B e ff/F e ff) 
used to derive the T mb scale are 0.45, 0.48 and 0.75 for the 12 CO 
1-0, 13 CO 2-1 and C' 8 Q 1- lines respectively faever et alJ 
1998l lKramer & Wild Il 994). All datasets have been smoothed 
to a resolution of 1' for comparison with our observations. 
These CO datasets provide deeper insight into the physical con- 

1 http://www.iram.fr/IRAMFR/GILDAS 
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Fig. 1. Velocity-integrated intensity maps of [Ci] 3 Pi- 3 Po emission (color) overlayed with contours of (a) 12 CO 4-3, 13 CO 
2—1, (c) 12 CO 1-0 and (d) C ls O 1-0 intensities. All maps are integrated over the range Vlsr=-20 to -5 km s . Colorscale of the 
[C i] map and contour levels of all other line intensities are indicated at the top left corner of each panel. The red star denotes the 
position of the hot spot. All data are shown at a common resolution of 1'. 



ditions of the regions emitting [C i] and CO 4-3 emission. 
12 CO 1-0 being optically thick, is most likely to be arising from 
closer to the PDR surface. The 13 CO 2-1 and C ls O 1-0 on the 
other hand are most likely to be optically thin and hence arise 
from deeper inside the clumps, and hence are crucial for as- 
certaining the physical parameters (like the optical extinction) 
characterizing the emitting regions. 

2.3. HIRES processed IRAS maps 

Dust continuum emission accounts for about 98% of the cool- 
ing of the molecular clouds. Here we have used FIR continuum 
maps to estimate the intensity of the FUV radiation field crucial 
for the chemical and physic al structure of the PD Rs. We have 
obtained HIRES-processed lAumann et al.lfl990l) IRAS maps 
of Cepheus B at 60 and 100 fim from the Infrared Processing 



and Analysis Center, Caltech. These dust continuum maps have 
angular resolutions (1'5) comparable to our sub-millimeter 
datasets. 

3. Observational Results 

3.1. Integrated Intensity Maps 

Figure [Qa) shows the integrated intensity maps of [C i] and 
12 CO 4-3 emission from Cepheus B. The maps are centered 
at a = 22 /! 57 m 25 s .7 6 = 62°34'40'.'7 (J2000). The 12 CO 4-3 
emission peaks at the position of the hotspot, and drops off sig- 
nificantly towards the south-east, and does not show much sub- 
structures. The [C i] emission in contrast though rather weak at 
the position of the hotspot, has almost uniform intensity over 
the ridge extended along the east-west direction with an iso- 
lated peak to the south-west (below the hotspot) and an ex- 
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tended high intensity peak deep into the molecular cloud to the 
east of the map. The somewhat different distribution of 12 CO 4— 
3 and [C i] emission suggest that while the 12 CO 4-3 emission 
traces higher temperature and densities as in the hotspot, the 
[C i] emission being optically thin traces the embedded PDR 
surfaces of the molecular material with higher column den- 
sity better. The features seen in the 12 CO 4-3 map agree rea- 
sonably well with the features seen in the 12 CO 3-2 map by 
iBeuther et alJfeOOOl) . 

Figures (b), (c) and (d) show overlays of the observed 
[Ci] emission with the 13 CO 2-1, 12 CO 1-0 and C 18 1-0 
intensities respectively. The [Ci] and 13 CO 2-1 appear to be 
almost concomitant tracing similar features in the Cepheus B 
cloud. To the west of the map the 12 CO 1-0 emission is similar 
in parts with the [C i] emission. However the 12 CO 1-0 intensi- 
ties drop off more drastically towards the east, which is similar 
to the features seen in the I2 CO 4-3 map. This is probably due 
to the higher column densities to the east, as manifested by the 
peaks in the 13 CO and C I8 maps. The lower I2 CO 4-3 emis- 
sion to the east, stem from a combination of lower densities and 
kinetic temperatures. The C ls O 1-0 emission, tracing higher 
column densities, looks markedly different from the 12 CO 4— 
3 emission and almost similar to the [C i] emission except for 
the region close to the hotspot where, the C ls O emission is al- 
most absent. The faint 13 CO 2-1 emission and almost no C 18 
1-0 emission from th e hotspot is consistent with the results of 
iMinchin et ail i 199 A who found a distinct lack of 13 CO emis- 
sion at the position of the hotspot. 

3.2. Velocity Structure 

Previous low-J CO observations showed that the Cepheus B 
emission consists of two major velocity components: one, the 
north-western part of the cloud closer to the hot spot emits at 

15 km s _1 , confirming the association of the hot spot with 

S155 having the same system velocity an d second, the south- 
eastern part emittin g at ~ -12 km s 1 (Beu ther et alJ feoOO: 
iMinchin et all 19 92). 

An overlay of the [C i] and 12 CO 4-3 channel maps (Fig.|2j 
shows that bulk of the cloud emits between -14 and -11 km s _1 , 
although substructures contributing to the 12 CO 4-3 emission 
change from the hotspot to the clump to the south-west of the 
map. The emission from the hotspot is slightly blue-shifted, be- 
tween -16 and -14 km s _1 . At -15.5 km s _1 , we detect in [Ci] 
a clump to the extreme north of the map. At -10.5 km s _1 the 
12 CO 4-3 emission from the cloud is too low and the [C i] emis- 
sion is restricted mainly to the east of the ridge. It is interesting 
to note that at -11.5 kms 1 the [Ci] and 12 CO 4-3 emission 
from the south-western clump are exactly concomitant. 

Figure[3]shows the [Ci], CO and [Cn] 158 yum spectra ob- 
served at the positio n of th e hotspot. The [C n] spectrum was 
observed by iBoreiko et alJ (l990) using the Kuiper Airborne 
Observatory with a resolution of 43", similar to the resolution 
of the other spectra (60"). The velocity resolution is 0.8 km s . 
It is one of the very few heterodyne spectra of the [C n] line 
taken so far. These spectra show that at the position of the 
hotspot the PDR emission is dominated by the velocity com- 
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Fig. 2. Overlay of [Ci] (color) and l2 CO 4-3 (contours) veloc- 
ity channel maps of Cepheus B. The colorscale ranges between 
1 and 1 1.6 K km s and the 12 CO 4-3 contours are between 5 
and 120 K km s in steps of 15 K km s" 1 . 



Table 2. Velocity components derived from gaussian fitting of 
spectra at the position of the hotspot in Cepheus B. 



Line 


Vlsr 
(km s" 1 ) 


FWHM 

(kms- 1 ) 


[Ci] 


-14.1 


1.9 


12 CO 4-3 


-13.9 


2.8 


12 CO 1-0 


-13.6 


2.2 


13 CO 2-1 


-13.4 


2.1 


[Cn] Comp 1 


-13.6 


3.8 


[C n] Comp 2 


-16.5 


6.0 



ponent at Vlsr of -13.8 km s _1 and the lineshapes with the ex- 
ception of [C n] are all primarily gaussian. The [C n] line at the 
hotspot position app ears to be the broade st, with traces of an 
extended blue wing. IBoreiko et al.l |[l590) have explained the 
[C n] emission from the hotspot position to be from two ve- 
locity components, the main narrow component is centered at 
-13.6 km s with a width of 3.8 km s _1 , while the blue compo- 
nent is centered at -16.8 km s _1 and has a width of 6.0 km s _1 . 
Table|3presents results of fitting Gaussian profiles to the differ- 
ent spectral lines observed at the hotspot. Broadly we find that 
both [Ci] and CO emission arise around Vlsr of -13.8 km s _1 
and have linewidths ~ 2.2 kms -1 . We have reproduced the 



Mookerjea et aL: PDRs in Cepheus B 



5 



lJ 


i 1 i i i i 1 i i 

-, (-145"/174") - 

CO 4-3 

U CI (x4) 

CII 

CO 1-0 

13 CO 2-1 

V>kjn Jl\ IT ._ 





-30 -20 -10 

Velocity (kra/s) 



Fig. 3. Observed CI, 12 C04-3, 12 COl-0, 13 C02-1 & [Cn] 
spectra at the position of the molecular hotsp ot in Ceph eus B. 
The [Cn] spectrum is adopted from Boreik o et alJ 111990). The 
ratio of integrated intensities 12 CO4-3/ 12 COl-0 at this position 
is 1.9. All data are shown at a common resolution of 1', except 
for the [C n] data which has a resolution of 43". 

gaussian fits to the [C n] spectra obtained by iBoreiko et alJ 
(1990) using the same parameters. The broader blue compo- 
nent at -16.5 km s _1 does not appear to be associated with the 
emission from the PDR/molecular cloud at all, rather the ve- 
locity matches the Vt„sr shown b y Ha emission from the H n 
region S155 llBoreiko et al.l l990). Based on the gaussian fit to 
the [Cn] spectrum, we estimate ~ 77% of the measured [Cn] 
emission arises from the PDR and the remaining from the ion- 
ized gas. The estimated relative contribution of PDRs to [Cn] 
emiss ion from nearby galaxies using the [N nl 1 22 fim is ~ 70- 
85% jMalhotraetalMoOll IKramer et alJ E)05). roughly con- 
sistent with the above result for Cepheus B. 

4. LTE Analysis 

4.1. The CO(4-3)/CO(1-0) distribution 

Figuregla) shows a map of the ratio of the 12 CO4-3/ 12 COl-0 
intensities integrated over the velocity range between -20 and 
-5 km s for the region in Cepheus B. It shows a smooth gra- 
dient between ratios of 0.3 - 0.6 in the south-east and ratios of 
1.5 - 2.2 to the north-west of the mapped region. 

Ratios of 0.9+0.2, as seen towards the bulk of the molec- 
ular ridge, suggest that both transitions are thermalized and 
can be explained in terms of LTE and optically thick emission. 
Neglecting the Rayleigh-Jeans correction, the ratio is expected 
to be 1 .0 in this case. The calibration uncertainty of line ratios 
is ~ 20%. The 4-3/1-0 ratio is sensitive to excitation temper- 
atures only upto 50 K, beyond that the ratio remains almost 
constant with increasing T ex . 

While ratios less than 0.8 suggest sub-thermal excitation 
of 12 C04-3, ratios above 1.2 are inconsistent with above as- 
sumptions. At the position of the hotspot for which the spectra 



are shown in Fig.|3]the integrated intensity ratio 4-3/1-0 has a 
value of 1.9. Similar to this position, at all other positions with 
high 4-3/1-0 ratios, we do not see any evidence of self absorp- 
tion of either line by cold foreground gas. Within the frame- 
work of single T ex and LTE, differences in optical thickness and 
volume filling factors of the emitting regions of the two lines 
could still partially explain the observed high 4-3/1-0 ratios. 
However we consider strong outwardly increasing temperature 
gradients to be the more plausible explanation. Such a gradient 
results in an inconsistency in deriving a single T ex from this ra- 
tio, since the 4-3 transition will predominantly trace the warm, 
still dense, outer regions, while the 1-0 traces the somewhat in- 
ner, colder regions. Such high I2 CO ratios are typically found 
in the s urface regions of molecular clouds subject to strong UV 
fields JCastets et al J 1 1 990t IKramer et all 1 1 996t iBeuther et all 
l2000UKramer et all2004lT 

Excitation temperatures derived from the 4-3/1-0 line ra- 
tios vary smoothly between ~ 8 K in the south-east of the 
mapped region and more than 50 K in the north-west near the 
position of the hot core and at the western edge of the molecular 
cloud (Fig.01b)). In the following discussion, we have assumed 
a constant excitation temperature of 50 K for4-3/l-0> 0.9. 

4.2. CP /CO Abundances in Cepheus B 

We have estimated the C° and CO column densities based on 
LTE approximation. For these calculations we have used the 
observed [Ci] and C 18 1-0 line intensities, the T ex estimated 
from the 12 CO4-3/ 12 COl-0 ratios and a relative abundance 
ratio [CO]/[C I8 0] of 500 JWilson & Rooc3ll994l) . Further we 
have assumed both [C i] and C 18 1-0 to be optically thin. The 
more detailed PDR analysis presented in Sect.|6j strongly sup- 
ports this latter assumption. 

Assuming LTE N(CO) depends linearly on T ex above 20 K, 
while in contrast, N(C) is rather independent of the assumed 
T ex to within 20% between 20 and 150 K. However, in between 
10 and 20 K, N(C) drops by more than a factor of 2, while 
N(CO) remains constant within 30%. 

Figure |3 a) shows a scatterplot of the the C° column den- 
sities, N(C) versus the CO column densities, N(CO). CO col- 
umn densities vary by more than one order of magnitude be- 
tween 10 17 and 2 10 18 cm -2 , corresponding to optical extinc- 
tions between ~ 1 mag and 30 mag. Although the range of val- 
ues of N(C) at different positions appear to be similar to the 
range shown by N(CO), a histogram of N(C) shows that for 
most positions N(C) does not show much variation, staying al- 
most constant at a mean value of 2.0 10 17 cirT 2 with an rms of 
1.410 17 cnr 2 . 

Figure EJb) shows a map of the N(C)/N(CO) (C/CO here- 
after) overlayed with contours of integrated intensities of C I8 0. 
This plot shows that regions of higher C I8 1-0 intensities, 
presumably tracing high H2 column densities, show rather low 
values of C/CO, in contrast to the outskirts. The ridge and the 
south-western core show low values of C/CO (~ 0.1), while 
the ratio rises gradually to values up to 4 outside these regions 
where the CO emission is more diffuse. 
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Fig. 5. (a) N(C) vs N(CO) for independent positions (positions separated by 1', i.e., beamwidth) in the map. Typical errors of 
15% for N(C) and 10% for N(CO) are shown for one position. (b)The C°/CO abundance ratio in Cepheus B overlayed with the 
contours of integrated intensity of C 18 Ol-0. We have considered here only those positions for which we have 12 C04-3 data and 
for which the [C i] and C ls O intensities lie above 3 times their respective rms noise levels. 



The constancy of C column densities is also reflected in 
Fig-IHshowing a scatterplot of the C/CO abundance ratio versus 
H2 column densities or optical extinctions derived from N(CO) 
assuming a CO/H 2 abundance of 8 1(T 5 (Blake et al. 1987). 
The C/CO ratio is ~ 1 for 2 mag and drops to 0.1 for 20 mag. 
A linear fit to the data results in 

log [n(C°)/N(CO)] = (-1.01 + 0.03W(H 2 ) + (21.5 + 0.72)(1) 

with a correlation coefficient of -0.82. The slope is almost 
-1 as would be expected for nearly constant N(C). 

Figure|6]also shows the C/CO abundances derived in other 
Galactic star forming regions and diffuse clouds. These stud- 



ies cover a range of optical extinctions between 0.3 mag and 
~ 100 mag with C/CO abundances varying by more than two 
orders of magnitude between ~ 4 and ~ 10~ 2 , approximately 
following the slope we find in Cepheus B. 

Clearly, the [C 1] 1-0 line is not a straightforward tracer 
of optical extinctions or total H2 column densities and total 
masses. 

The scatter of C/CO abundance varies with Ay. For Ay's 
upto 20 mag the scatter decreases from being almost compa- 
rable to the mean C/CO itself to about 50% of the mean. At 
~ 20 mag the scatter is the smallest, with the C/CO abundance 
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from five different data sets agreeing to a value of 0.1 with 
a scatter of 0.05. Note that several studies which cover only 
a small range of extinctions find rather constant C/CO abun- 
dances. 

In Sect. [6] we show that clumpy PDR models naturally ex- 
plain the observed variation of the C/CO abundance. 
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Fig. 6. Column density ratio of C to CO in Cepheus B. 
Plotted as in Fig |5] are independent positions, which are 
separated by 1'. Also shown are observed valu es in ot her 
Galac tic mo lecular clouds like NGC 7 023 JGerin et alJ 
19981) . S106 dSchneider et al.ll2003h. PR 15 dOka et alJl200ll) . 
MCLP 123.5+24.9 JBensch et alJ20o"3h. High Latitud e Clouds 
i Ingalls et alJ 1 1 9971). NGC 1333 JOka eTail l2004. o Oph 
i Kamegai et alJ I2003I). W3Main jKramer et all l2004 and 
Orion Jlkeda et alJl2002l). The points corresponding to GMCs 
are taken from P lume et all i 19991) in which [C i] data of Cep A, 
S 140, NGC 2024 and W3 were presented. The red drawn 
line corresponds to the parameter derived by fitting only the 
Cepheus B data. Typical errors of 15% for N(C)/N(CO) and 
10% for N(H2) are shown for one position. The blue line is the 
result of fitting data points for all previous observations. 



5. FUV intensity from FIR continuum 

The spectral lines observed in Cepheus B which are used here 
are not sensitive to the FUV radiation field, and there are no 
largescale [C n] observations available in the literature for this 
region. We thus estimate the FUV radiation field from the ob- 
served FIR continuum. Because dust grains absorb most of the 
FUV photons incident on the cloud and most of that energy 
is reemitted in the form of FIR continuum, and other heating 
mechanisms like cosmic ray heating can be neglected, the in- 
tensity of the FUV radiation can be inferred from the FIR con- 
tinuum. This estimate is valid even if a fraction of FUV photons 
escape without impinging the cloud. 



We have estimated the FIR intensity of the mapped region 
using the HIRES proce ssed 60 and 100 um I RAS fluxes fol- 
lowing the relationship (Nakagaw a et alll998l) : 



I FIR =3.25x10" 



/ v (60^m) 



Jy sr 



+ 1.26x10" 



/ v (100//m) 



Jy sr 



(2) 



The FUV intensity expressed in u nits of the Praine field 
dPrainel 1 978t IPraine & Bertoldi 1996) is given by 



1 



An 



X 1.71(2x1.610- 3 ) /fiR (3) 

Following Kaufm an et alJ Jl999l) . this formula takes into 
account that there is about equal heating of the grains by pho- 
tons outside the FUV band (leading to the factor 2) and that the 
bolometric dust continuum emission is a factor of ~ 2 larger 
than the FIR intensity derived above from the IRAS fluxes. We 
also assume that the emission stems from slabs which are il- 
luminated from all sides and the FIR emission stems from the 
back and front sides. The latter two corrections approximately 
cancel out. Finally, in order to convert the FUV field e xpressed 
in un its of Habing field (1.6xl0" 3 ergs" 1 cm" 2 sr" 1 ; Habing 
1968) to Praine field it is required to divide by 1.71. 

(arcsec) 



Contours X (in Draine Unit) 
25 50 75 100 200 300 400 600 800 1000 1400 G 




a (2000) 

Fig. 7. FUV intensity distribution (red contours) in Cepheus B 
estimated from the FIR fluxes measured by IRAS, overlayed 
with [C i] integrated intensities (in color). Contour levels of the 
FUV intensity are indicated at the top left corner. The intensity 
scale for the [Ci] map are the same as in Fig.^ White stars 
denote HP2 17086 (north) and HP2 17061, Crosses denote the 
NIR sources A-NIR (south) and B-NIR and the filled triangle 
denotes the location of the hot spot. 



FigureQshows the distribution of FUV intensity calculated 
as above, overlayed with the observed [Ci] emission. Within 
the region mapped, the FUV field varies smoothly between ~ 
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25 in the south-east and 1500 times Draine field, in the north- 
west near the hot spot. This variation of x correlates well with 
the corresponding variation of excitation temperatures derived 
from the 12 CO line ratios (Fig.|3b)). 

In Fig. we have marked the positions of the nearest O 
and B stars, HD2 17086 (north-west), H D2 17061 (west) an d 
the embedded NIR sources A and B (cf. iTesti etaD (.1995)). 
NIR-A, which does not have any infrared excess and which is 
deemed to be a highly reddened main sequence star of spectral 
type B 1-B0.5, is responsible for the excitation of the associated 
blister Hn region (the radio continuum source A). The FUV 
intensity distribution calculated as above is consistent with the 
results of Test i et al.l II 1995ft and suggest that the star A-NIR is 
mostly responsible for the FIR emission from the region. The 
total estimated FIR luminosity from the IRAS fluxes and for 
the assumed distance of 730 pc for Cepheus B is ~ 10 4 L and 
this is consistent with the luminosity due to a Bl type star as 
well. 

6. PDR Model interpretation of C°/CO vs N(H 2 ) 
relation 

We interpret here the observed anti-correlation of the C/CO 
abundance ratio with the hydrogen column density (N(H2)) in 
terms of PDRs created due to the FUV irradiation of molecu- 
lar material. For this we have used the K OSMA-r PDR model 

i h 994 r 



^rzer et alJ 1199 6: Rollig et alJ 120061) . This model consid- 
ers spherical clouds illuminated by an isotropic FUV radia- 
tion field and cosmic rays, calculates the resulting chemical 
and temperature structure of the cloud. This is followed by 
radiation transfer calculations to estimate the emission of dif- 
ferent species. The KOSMA-r model has 13 CO chemistry in- 
cluded in its network. The emission from the models are cal- 
culated as a function of the hydrogen volume density, FUV ra- 
diation field and mass of the spherical clump. The KOSMA-r 
model assumes that a clump has a power-law density profile of 
n(r) ~ r~ L5 for0.2< r/r c / < 1 and «(r) = const, for rjr c \ < 0.2. 
The surface density n c \ is about half the average clump density. 

6.1. Estimates of volume density 

First, we have used PDR models to derive the range of allowed 
hydrogen volume densities based on the observed CI/ 13 CO 2-1 
and CI/ 12 C04-3 intensity ratios. At all the positions for which 
we have calculated C/CO abundance ratio the CI/ 13 CO 2-1 
ratio varies between 0.3 and 0.7, while the CI/ 12 C04-3 ratio 
varies between 0.2 and 0.6. Based on FIR maps the FUV field 
in the mapped region varies between 50 and 10 3 Draine units. 
Figure [8] shows the contours bracketing the observed two line 
ratios for constant FUV field and different PDR clump densi- 
ties and masses. These two line ratios do not allow to discrim- 
inate between different clump masses. The best fitting clump 
densities lie between 2 10 4 and 8 10 5 cm" 3 , for clump masses 
between 10~ 2 and 10 M . This range of densities is rather in- 
dependent of the FUV field. We note that the clump masses are 
not very well constrained, however the intersection of the two 
ratio contours in Fig.[S]provide physically reasonable limits on 
clump masses, which have been loosely used in our subsequent 
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Fig. 8. Results from the PDR models for external FUV fields 
of x - 10 3 (a) and x = 10 2 , expressed in Draine units (b). 
Contours show integrated intensity line ratios of CI/ 12 C04-3 
and CI/ I3 C02-1. The range of observed values is shown: 0.3- 
0.7 for CI/ I3 CO 2-1, 0.2-0.6 for CI/ 12 C04-3. We consider 
only positions for which the C/CO abundance ratio has been 
estimated. Both ratios suggest that the volume densities on an 
average lie between 10 4 and 10 5 cm 4 . 



analysis, but do not have any significant impact on our final re- 
sults. The observed I2 CO 4-3/1-0 line ratios are consistent with 
this interpretation. 

6.2. Discussion 

6.2.1. Single Spherical PDR Model 



10 



10 



c 



10 



10 



-i 1 1 1- 



Mci U m P =lM« n Rtot -0.037pc 
«clump=0 001«»» Riot = 0.0037pc 




CO 



n=10 cm X=l° 



0.5 

Relative Depth d/R lot 



Fig. 9. Abundance profiles of C° and CO calculated by the 
KOSMA-r PDR model for two clumps with masses 1 M and 
0.001 M . The x-axis is given in units of the relative depth 
d/R lot . A value of 1 denotes the center of the cloud. The FUV 
field strength is 10 2 in Draine units. 



We first explore the single spherical PDR models described 
above, in order to derive the observed C/CO vs N(H2) anti- 
correlation. As shown above, clumps with masses between 
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10 2 and 10 M with densities ~ 10 5 cm 3 reproduce the range 
of [C i]/ 13 C02-l and [C i]/ 12 C04-3 intensity ratios observed in 
Cepheus B. In Fig.|5J we present the calculated radial variation 
of the C and CO abundances for two PDR clumps with masses 
of 10~ 3 and 1 M . The surface density is fixed at 10 5 cm~ 3 and 
the FUV field is fixed at 10 2 (in Draine units), the average value 
derived from the HIRES data in Sect. [5] 

At the surface of the PDRs carbon is fully ionized by the 
incident FUV radiation. Increasing depth leads to increased 
absorption of the FUV radiation by both gas and dust. This 
allows C + to recombine with electrons to form C°. At even 
larger depths CO starts to form because the photodissociating 
FUV photons are attenuated sufficiently leading to a chemi- 
cal stratification typical of PDR models, C + /C°/CO The ab- 
solute depth where this transition occurs depends on the gas 
density n and the FUV intensity x- In spherical models it also 
depends on the size R tot of the cloud since the local mean inten- 
sity J v (r) = (AnY 1 J I v (r, Q) dQ, at a certain depth is larger for 
a small clump compared to a bigger one. Accordingly a clump 
with M c i = 1 Mq consists of a much larger CO core than a 
clump with M c i = 0.001 M Q . The average CO column density 
increases with mass. The dependence of N(C) on the mass is 
weaker. Once the clump becomes massive enough to signifi- 
cantly shield the photodissociating FUV photons a further in- 
crease in mass does not lead to larger C° columns but larger CO 
column densities. This behavior is illustrated in Fig. [9] with the 
abundance of CO at the inner core (with relative depth=l) of 
the clumps being identical in both cases. Thus the 1 M clump 
has a much higher total fraction of CO and lower C/CO than the 
lower mass clump. This indicates that the lines of sight hav- 
ing high C/CO abundance ratios are dominated by clumps of 
smaller sizes. 

Table 3. KOSMA-r calculations of H 2 , C and CO column 
densities of spherical clumps with n=10 5 cm' 3 and^f = 10 2 
(Draine units). The column densities are for individual clumps, 
and not the column densities observed per beam. 
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6.6 


14 
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0.004 



Table[3]lists the clump-averaged C, CO and H2 column den- 
sities of the KOSMA-t PDR model, considering a single spher- 
ical clump. We note here that the column densities in Table 
correspond to those calculated by the PDR model for a single 
clump and do not corroborate with the observed column densi- 
ties per beam, as plotted in Fig. [6] The C/CO ratios calculated 
by the PDR model is correlated with the mass of the clump. We 
have checked the peak optical depths of the [C 1] line emission 
in a pencil beam centered on the clump and for the line cen- 



ter position. For all clumps studied here, these peak opacities 
are less than 0.5. The clump and line averaged opacities of the 
[C 1] line are much less. This justifies our previous assumption, 
used in the LTE analysis, that the [Ci] line is optically thin. 
Similarly, the modelled 13 CO peak opacities indicate that the 
C ls O line is optically thin. 

PDR models with n=10 5 cirT 3 have clump sizes ranging 
between 0.004 pc (for m c i ump = 10~ 3 M ) and 0.08 pc (for 
rriciump = 10 M ). In comparison, a beam of 60" corresponds 
to a diameter of 0.21 pc at the distance of Cepheus B. It is thus 
not possible to compare the observed C/CO vs N(H2) per 60" 
beam correlation directly with the N(H2) per clump (N(H2) c i in 
Tabled} obtained from the single spherical PDR clumps. 

6.2.2. The Clumpy PDR Scenario 

The observed uniformity of the [Ci] intensity distribution 
throughout the PDR (Fig. and the failure to interpret the 
C/CO vs N(H2) using single spherical PDR model leaves open 
two possible physical scenarios: (1) dumpiness of the medium, 
i.e., either there are multiple plane-parallel slabs stacked or 
there is an ensemble of clumps along the line of sight, and 
(ii) a single, homogeneous plane-parallel PDR slab inclined 
relative to the major sources of UV radiation HD217086 and 
HD217061. 



Mass in 1' beam (M^J^j 




N(H a ) [em" 8 ] 



Fig. 10. KOSMA-r PDR model interpretation of the observed 
C/CO vs N(H2) anti-correlation. For comparison results of 
face-on homogeneous plane-parallel models are also shown. 
Drawn lines correspond to models the input parameters for 
which are marked at the bottom left corner. M c i corresponds to 
mass of an individual clump. Mb eam shown along the upper x- 
axis is the mass of gas in the 60" beam. The number of clumps 
in the beam is given by Mb ea m/M c i. 
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We first consider the simplified clumpy PDR interpretation 
using ensemble of clumps having the same mass and density at 
each position. Figure ITU1 shows that the observed C/CO abun- 
dances and total H2 column densities of Cepheus B can be 
reproduced with clumpy PDR models. It shows results based 
on the ensemble interpretation using the spherical PDR mod- 
els with n=10 5 cm 3 , x — 10 3 and 10 2 and n=10 4 cirT 3 and 
X - 10 3 . The drawn lines correspond to model calculations. We 
consider an ensemble of clumps with mass M c i (also marked 
in Fig. II Oi contributing to a total mass of Mt eam in a 60" 
beam. The total number of clumps in the beam is given by 
N en s= Mbeam/Mci- Table |4] summarizes the defining parame- 
ters of the clump ensembles for which the C/CO abundance 
ratios are shown in Fig. [TO] We find that the observed C/CO 
anti-correlation depends primarily on the characteristic mass 
of individual clump along a particular line of sight and does 
not have a strong dependence on the assumed FUV field. For 
clumps with densities of 10 4 cm -3 , masses below 0. 1 M result 
in much higher values of C/CO ratios than observed. The area 
filling factors (</>a) of the ensemble vary between 1.6 and 3.3 
(Table0}, from the higher to the lower values of clump masses 
considered in our calculations. This indicates that the clumps 
fill the beam and in most cases partially overlap. 



Table 4. Results of PDR model calculations for ensembles of clumps 
having same mass and density. 0a is the area filling factor of the 
clumps. 
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* In Draine units 



For the sake of completeness we also explored the sec- 
ond possible physical scenario, an inclined plane-parallel PDR, 
which might produce the observed uniform PDR emission 
apparently deep into the cloud and the observed C/CO and 
N(H2) anti-correlation. For this purpose we have used the 
plan e-parallel PDR models illuminated from one side by 
Koest er et al .I d 1994 in order to explicitly calculate the C/CO 
ratios and N(H2) for selected values of volume densities and 
UV-field. We have explored several values of «, slab-size 
(N(H2)) and x an d nave found consistently that for face-on 
PDRs the calculated C/CO ratio drops much faster with N(H2) 



as compared to the observations. Figure^|shows the results of 
calculation for a plane-parallel model with n = 2.5 10 4 cm -3 
and^ = 100, for different slab-sizes or total NXH2). Changing 
the inclination of the PDR has the same effect on both N(C) 
and N(CO), and thus the C/CO ratio remains unaltered, though 
N(H2) changes as the reciprocal of the sinus of the inclina- 
tion angle. For this particular set of models we find that at 
N(H 2 ) = 3 10 21 cm -2 the calculated C/CO ratio is 0.08. In 
contrast, observationally we find the C/CO ratio to be 0.08 at 
NCH2) ~ 3 10 22 crrT 2 . There is a discrepancy by a factor of 10 
in N(H2), which can be explained by an inclination angle of 
6°. Similarly, different inclination angles are required to bring 
the different plane-parallel models to conform to the observed 
slope of the anti-correlation. Hence, though it might be possi- 
ble to explain individual positions by making suitable choice of 
inclination angles, it is not possible to reproduce the observed 
anti-correlation at all positions in Cepheus B using a single in- 
clination angle. Further, the anti-correlation as has been shown 
in Fig. [6] holds for a wide range of sources so it is rather un- 
likely that all of them would either have the same inclination 
angle or have different inclination angles and still conform to 
the same anti-correlation. The plane-parallel models thus fail to 
reproduce the slope of the C/CO vs NCH2) anti-correlation en- 
tirely which extends over two orders of magnitude in NXH2). 
Intuitively, both uniformity of [Ci] emission over the entire 
ridge of Cepheus B and the C/CO vs N(H2) anti-correlation 
seen in most PDRs appear to be obvious outcomes of homo- 
geneous plane-parallel models which are inclined. However, 
explicit calculations suggest that although the homogeneous 
plane parallel models predict almost constant N(C) with chang- 
ing NXH2), they do not reproduce the observed global C/CO vs 
N(H2) anti-correlation even when inclined. 

We conclude that the spherical clumpy PDR model repro- 
duces the frequently observed global C/CO vs N(H2) anti- 
correlation over the entire range of observed NCH2). These 
clumpy PDR models can thus be used to understand the typ- 
ical density and size (hence mass) of the regions that domi- 
nate the observed emission. This does not imply the absence 
of clumps of sizes smaller or larger than those reproducing the 
observed C/CO ratios. However the emission cannot be domi- 
nated by such clumps if the intensity and abundance ratios of 
these clumps differ significantly from the observed numbers. 
Hence the low observed relative abundance ratio C/CO origi- 
nate from positions of high column density where the emission 
is dominated by a few more massive clumps, while high ratios 
originate at positions covered by many small clumps, having a 
low line of sight column density. 

The calculations presented here are meant to be of qualita- 
tive nature in order to derive insight into the clump parameters 
which determine the observed C/CO ratio. 

Figure [5] shows that the C/CO ratio is higher towards the 
south-eastern part of the ridge further away from the sources 
of FUV radiation. This together with the results of clumpy 
PDR models which find that higher C/CO ratios are associated 
with clumps of smaller mass present an interesting situation. 
In Cepheus B towards the north-western edge there are evi- 
dences of embedded sources while there are none to the south- 
east. Thus, using the clumpy PDR models we see that regions 
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with embedded sources correspond to lines of sight with higher 
mass clumps. 



7. Summary 

We have presented the first largescale fully-sampled maps 
of [Ci] (492 GHz) and 12 CO 4-3, of the edge-on PDR, 
Cepheus B. The [Ci] emission does not delineate the PDR 
clearly, rather it traces the bulk emission of the more quies- 
cent part of the molecular ridge extending to the east. Bulk of 
the PDR gas emits at velocities between -14 and -11 km s _1 . 
The 12 CO 4-3 emission is strongly dominated by the hotspot 
present toward the north-western tip of the molecular cloud. 
We observed 12 CO4-3/l-0 intensity ratios in excess of 1.5 to- 
wards the north-western edge of the cloud. Such ratios can only 
be explained in terms of temperature gradients along the line of 
sight. Excitation temperature, derived from the observed 12 CO 
4-3/1-0 ratio, varies between 20 K (in the quiescent parts of 
the cloud) to 50 K (close to the hotspot). We conclude that [C i] 
primarily traces the higher column density regions while 12 CO 
4-3 traces the higher temperature and density regions. 

In Cepheus B we find that N(C) remains nearly constant 
although N(CO) varies over one order of magnitude. The ob- 
served C/CO ratio varies between 0.06 and 4.0 and shows a di- 
rect anti-correlation with the hydrogen column density N(H2). 

As estimated from the FIR continuum emission, the FUV 
radiation field in Cepheus B varies between 25 & 1500 in 
Draine units, with an average value of 100 over most of the 
[C i] emitting region. We have used the KOSMA-t PDR model 
to explain the observed intensity ratios and the observed C/CO 
vs N(H 2 ) anti-correlation. These PDR model calculations show 
that for^ = 100, the observed CI/ 13 C02-1 and CI/ 12 C04-3 
ratios indicate local hydrogen densities between 10 4 -10 5 cirT 3 
with even higher values at a few positions. The 12 CO 4-3/1-0 
ratio can be reproduced by the PDR models and the modelled 
clumps indeed exhibit outwardly increasing temperature gradi- 
ents, as predicted from LTE analysis. 

The observed global C/CO vs N(Ha) anti-correlation and 
its slope cannot be reproduced by single spherical PDR models 
or by plane-parallel homogeneous slabs with or without incli- 
nation. We were able to explain the observed anti-correlation 
between C/CO and N(H 2 ) considering an ensemble of clumps 
of the same mass and density (n=10 5 crrT 3 ) for each observed 
position. We find that (i) the C/CO ratio is characteristic of the 
dominant mass in the ensemble, clumps of smaller mass have 
higher C/CO ratio and (ii) clumps of larger mass have to neces- 
sarily dominate emission along lines of sight with higher col- 
umn density in order to reproduce the observed lower C/CO 
ratio. 
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